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Abstract —This paper investigates simultaneous wireless infor¬ 
mation and power transfer (SWIPT) for multiuser multiple- 
input-single-output (MISO) secrecy channel. First, transmit 
beamfoming without artificial noise (AN) design is considered, 
where two secrecy rate optimization frameworks (i.e., secrecy 
rate maximization and harvested energy maximization) are in¬ 
vestigated. These two optimization problems are not convex, and 
cannot be solved directly. For secrecy rate maximization problem, 
we employ bisection method to optimize the associated power 
minimization problem, and first-order Taylor series expansion is 
consider to approximate the energy harvesting (EH) constraint 
and the harvested energy maximization problem. Moreover, we 
extend our proposed algorithm to the associated robust schemes 
by incorporating with channel uncertainties, where two-level 
method is proposed for the harvested energy maximization 
problem. Then, transmit beamforming with AN design is studied 
for the same secrecy rate maximization problem, which are 
reformulated into semidefinite programming (SDP) based on one¬ 
dimensional search and successive convex approximation (SCA), 
respectively. Moreover, tightness analysis of rank relaxation is 
provided to show the optimal transmit covariance matrix exactly 
returns rank-one. Simulation results is provided to validate the 
performance of the proposed algorithm. 

Index Terms —MISO system, SWIPT, physical-layer secrecy, 
bisection method, successive convex approximation (SCA). 

I. Introduction 

As one of the main techniques in fifth-generation (5G) 
wireless networks, energy harvesting is a technique to extract 
energy from the external natural environment (i.e., solar power, 
wind energy, etc.), which can be applied to circumvent energy- 
limited issue and improve the energy efficiency of wireless 
networks m, 13 . However, traditional energy harvesting tech¬ 
niques are dependent on external sources, which are not part 
of the communication network. Since it is not always possible 
that the nature environment can provide a stable energy to the 
wireless devices, the energy can be alternatively provided by 
wireless devices, since radio frequency (RF) signals can carry 
energy that is used as a vehicle for transporting information 
0. In particular, the transmitter in a communication network 
not only can transmit the signal, but also transfer the power 
to the receivers who are employed to harvest energy to charge 
for their devices. Wireless power transfer has been a promis¬ 
ing paradigm to provide power supplies for communication 
devices to mitigate the energy scarcity and extend the lifetime 
of wireless networks 0, EL 

Secrecy performance played more and more important role 
in wireless networks, and there is an increasing number of 
interests focusing on physical-layer security. Unlike traditional 
cryptographic methods, which is employed to improve secrecy 
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performance in the network layer, physical-layer security was 
developed in terms of information-theoretical aspects improve 
the security of wireless transmission 0-0. Multiple an¬ 
tenna wiretap channel has been widely considered with the 
advantage of having extra degrees of freedom and diversity 
gains 0- In addition, some state-of-art techniques have been 
developed for multiple-antenna transceivers, which aims to 
introduce more interference to eavesdroppers, including coop¬ 
erative beamforming (CB), cooperative jammer (CJ), artificial 
noise (AN), energy harvesting beamforming (EHB), device- 
to-device transmission (D2D), etc q. Ha-in. The convex 
optimization techniques have been employed to solve the 
secrecy rate optimization problem (i.e., power minimization 
and secrecy rate maximization) based on perfect and imper¬ 
fect channel state information (CSI) IfTOj . The cooperative 
beamforming is employed that requires the relays to forward 
the signal from the source to the legitimate user based on 
the assumption that the direct transmission is not available. 
The optimal power allocation for single-relay multicarrier case 
based on decode-and-forward (DF) has been proposed to maxi¬ 
mize the sum secrecy rate El, whereas in E3 the relay works 
in amplify-and-forward (AF) scheme for MIMO system, where 
the source and relay beamforming is jointly designed to max¬ 
imize the secrecy rate in the cooperative scheme. Cooperative 
jamming is applied in improving physical layer security E2- 
el For single-antenna case, the secrecy rate is maximized 
by employing one-dimensional search algorithm El. In Oil, 
Taylor series expansion is investigated to approximate the 
secrecy rate optimization for MIMO channel with cooperative 
jammer, which can be reformulate the non-convex secrecy 
rate optimization framework to convex form, in addition, 
game theory technique (i.e., Stackelberg game) is applied for 
improving secrecy rate, and Stackelberg game equilibrium 
can be derived in terms of closed-form solutions. For MIMO 
relay networks, the optimal power allocation is derived for 
the generalized singular value decomposition (GSVD)-based 
secure relaying schemes G3L Joint optimization of transmit 
beamforming with artificial noise (AN) is proposed in |fl6l , 
whereas the energy harvesting beamforming (EHB) and D2D 
transmission are not only improve the secrecy rate, but also 
can harvest the energy to prolong the lifetime of the devices 
by recharging the batteries and enhancing D2D transmission 
rate, respectively E3-IH. 

More and more recent attentions are focusing on the com¬ 
bination of transmit beamforing with AN, where AN can be 
employed to mask the transmit beamforming so that improve 
the secrecy rate, in addition, the transmit beamforming and 
AN can be harvested by the energy harvesting (EH) receivers 
to improve energy harvesting performance on, Gl, El- 
Moreover, the eavesdroppers in these two literatures play the 


two roles: a) passive eavesdropping, b) energy harvesting. 
The secrecy rate maximization problem are formulated into 
semidefinite programming (SDP) and provide the proof of the 
rank-one for the transmit covariance matrix. Robust secure 
transmission for MISO Simultaneous Wireless Information 
and Power Transfer (SWIPT) system are proposed without 
AN ED and with AN E2, respectively by incorporating 
with channel uncertainties, where the proofs of the rank-two 
solution for the transmit covariance matrix have been provided. 
Motivated by the existing literatures ED, E2, this paper 
considers SWIPT for MISO secrecy system with multiple 
eavesdroppers and multiple EH receivers, here the secure 
transmit beamforming without AN and with AN are designed. 
The contribution of this paper is that we address the following 
problems: 

1) Secrecy rate optimization without AN: We formulate 
the secrecy rate optimization framework to design the 
transmit beamforming without AN. Two optimization 
problems (i.e., secrecy rate maximization and harvested 
energy maximization) are consider, which are not convex 
and cannot be solved directly. In order to circumvent 
this non-convexity issue, we propose a novel refor¬ 
mulation for the secrecy rate maximization problem 
based on bisection method, and the harvested energy 
maximization problem via first-order Taylor series based 
one-dimensional line search algorithm, respectively. In 
addition, these associated robust schemes are solved 
by incorporating with channel uncertainties. Unlike the 
most existing literature, which employs semidefinite 
relaxation (SDR) (or rank-relaxation) to design the trans¬ 
mit beamforming, in this paper, we consider the relax¬ 
ation to circumvent the rank-relaxation and we show 
our proposed schemes can achieve the same secrecy 
performance with that via SDR. 

2) Secrecy rate optimization with AN: Secrecy rate max¬ 
imization problem is formulated to jointly optimiza¬ 
tion for transmit beamforming with AN. In this paper, 
we propose two reformulations for this optimization 
framework based on one-dimensional line search and 
SCA algorithms, respectively, which can be used to 
reformulate the problem to semidefinite programming 
(SDP). In addition, we provide the tightness analysis for 
the rank-relaxation, which shows the optimal solution of 
the secrecy rate maximization exactly returns rank-one. 
Moreover, we extend our proposed schemes to the asso¬ 
ciated case by incorporating with channel uncertainties 
and the case with multiantenna eavesdroppers and EH 
receivers. 

The remaining part of this paper is organized as follows. A 
system model is described in Section [ID The secrecy rate opti¬ 
mization problems with SWIPT are solved to jointly optimize 
the transmit beaforming without AN based on perfect and 
imperfect CSI in Section [In] and [IV] Secrecy rate optimization 
problem with SWIPT to joint optimization of the transmit 
beaforming with AN is investigated in Section [V] Section 
[Vll provides the simulation results to support the proposed 
schemes. Finally, conclusions are presented in Section [VTIl 


A. Notation 

We use the upper case boldface letters for matrices and 
lower case boldface letters for vectors. (-) T , (•)* and (■) H 
denote the transpose, conjugate and conjugate transpose re¬ 
spectively. Tr(-) and E{-} stand for trace of a matrix and 
the statistical expectation for random variables. Vec(A) is the 
vector obtained by stacking the columns of A on top of one 
another and (xi is the Kronecker product. A A 0 indicates 
that A is a positive semidefinite matrix. I and (-) -1 denote 
the identity matrix with appropriate size and the inverse of a 
matrix respectively. || • ||2 represents the Euclidean norm of 
a matrix. ‘ft{ ■} stands for the real part of a complex number, 
whereas |A| denotes the determinant of A. [a;] + represents 
max{a;, 0}. 


II. System Model 

We consider a MISO secured SWIPT channel, where it con¬ 
sists of one multiantenna legitimate transmitter, one legitimate 
user, K eavesdroppers and L energy harvesting (EH) receiver. 
It is assumed that the transmitter is equipped with Nt transmit 
antennas, whereas the legitimate user, the eavesdroppers and 
the EH receivers each have a single receive antenna. The 
channel coefficients between the legitimate transmitter and the 
legitimate user, the k -th eavesdropper as well as the /-th EH 
receiver are denoted by h s £ C^ TXl , h e ^ £ C Ntx1 , and 
h£ C NrXl , respectively. The noise power at the legitimate 
user and the eavesdroppers are assumed to be rA and er 2 . The 
received signal at the legitimate user and the fc-thl eavesdropper 
can be written as 


y s = hf ws + n s , 

y ek = hej fc ws+n ej fc, k=l,...,K, 

where s(E{s 2 } = 1) and w £ <C Ntx1 are the desired signal 
for the legitimate user and the beamforming at the legitimate 
transmitter, respectively. In addition, n, ~ CA'"(0, cr 2 ) and 
n e ,fc ~ C,V"(0, ct 2 ) represent the noise of the legitimate user 
and the fc-th eavesdropper, respectively. Thus, the achieved 
secrecy rate at the legitimate user is expressed as follows: 




max log 
k 



, Vfc. 


(i) 


The harvested energy at the /-the EH receiver is written as 
Et = |hfw| 2 , V/. (2) 


Remark: This system model is considered to consists of L EH 
receivers, which harvest a number amount of power carried 
by RF signal without AN or with AN based on a reliable 
transmission scenario. These EH receivers sometimes play a “ 
helper ” role by employing the harvested power to introduce a 
jamming signal to confuse the eavesdroppers ll23l . However, 
the efficiency of this harvest-and-jamming policy is dependant 
on the network topology |24l . In this paper, we mainly focus 
on the secrecy rate optimization without AN or with AN, 
based on this reliable communication, the performance the EH 
receiver harvest the power is exploited, whereas how to use 
the harvested power is beyond our paper. 






TABLE I: Bisection methods 

1) Given lower and upper bound of the targeted secrecy rate R m - in and 
i?max, and a desired solution accuracy r (very small value). 

2) Setting R = (R m in + -Rmax)/2. 

3) Iteration loop begin 

a) Solve the corresponding power minimization problem in s 
using the relaxation method to obtain the beamforming w. 

b) Computing the transmit power P = ||w|| 2 . 

c) If P < P, then Rmin = R', otherwise, Rm ax = R. 

d) Until R max /i In i T1 < r, break. 

4) Iteration loop end 

5) R is the achieved secrecy rate of the secrecy rate maximization 
problem, and w is the corresponding optimal solution. 


III. Secrecy Rate Optimization without AN 
In this section, we consider the MISO secrecy with multiple 
eavesdroppers and EH receivers shown in Section QI] where 
the secure beamforming is designed without AN, in addition, 
the secrecy rate optimization frameworks (i.e., secrecy rate 
maximization and harvested energy maximization) are refor¬ 
mulated as a novel SDP form to circumvent rank relaxation. 
Thus, these optimization problems are formulated as follows: 

• Secrecy rate maximization: The secrecy rate is maximized 
subject to transmit power and minimum harvested energy 
constraints. 

maxminP s , s.t. ||w11“ < P, mini?; > E, Vk,l, (3) 

w k l 

where P is the maximum available transmit power at 
the legitimate transmitter, E denotes the target harvested 
energy of the EH receivers. 

• Harvested energy maximization: The harvested energy is 
maximized with the constraints of the secrecy rate and 
transmit power. 

maxminSi, s.t. minii 3i j. > R, ||w|| 2 < P, Vk,l, (4) 

w l k ’ 

where R is the predefined secrecy rate of the legitimate 
user. 

These above problems are not convex, and cannot be solved 
directly. Unlike existing literature EH . where rank relaxation 
is consider to reformulate the secrecy rate, however, since it 
is challenging to show rank-one solution, the authors provided 
suboptimal solution for the secrecy rate maximization. In this 
paper, we propose a novel relax method for secrecy rate 
constraint to circumvent rank relaxation. 

A. Secrecy Rate Maximization 

In this subsection, we consider the relax of the secrecy rate 
maximization problem £3}, which is not convex due to the se¬ 
crecy rate and non-convex EH constraint, and cannot be solved 
directly. In order to circumvent these two issues, we convert 
this problem into a sequence of the power minimization 
problems, one for each target rate R > 0. The optimal solution 
of the secrecy rate maximization problem can be obtained by 
solving the corresponding power minimization problem with 
different R, which can be obtained by using bisection search 
over R. Thus, the bisection method can be summarized in 
Tab I eQ] to solve this secrecy rate maximization problem. In the 
following, we will focus on the power minimization problem 
written as follows: 

II ll 2 

mm ||w|| , 

W 

s.t. minP s > R, min Ei > E, Vk,l. (5) 

k l 


Now, we consider the power minimization problem based 
on the assumption that the transmitter has perfect CSI of the 
legitimate user, the eavesdroppers and the EH receivers. Thus, 
the problem in 0 can be relaxed as 


min 

W 

S.t. 


log 1 + 


-log 14- 


|h, w f>E, VI. 


( 6 ) 


The above problem is not convex in terms of w and the 
non-convex EH constraint. In order to circumvent these issues, 
we convert the secrecy rate constraint into SDP form and 
reformulate the EH constraint based on successive convex 
approximation (SCA), respectively. Thus, the following 
lemma is required to tackle the power minimization problem 
©: 


Lemma 1: The problem in 0 is reformulated into the 
following form: 


mm si 

W 

S.t. 


hx o, 


St = 


w ff h s I 


R 
2T 


w h e ,fc 

R 


^-w ff h e , fc 
e i 
(2 R —1)2 


-w h s 


y o, Vk, 


(2 k -1)4 

xi = 5R{w w h;}, yi = 9{w ff h;}, u; = [xi yi\ 


,(")||2 _ 


+ 2 J2 u i (* *) \. u ‘ (*) - u i (*)] > vz. 


(7) 


Proof: Please refer to Appendix [A] ■ 

In problem ([7}, the secrecy rate constraint is reformulated 
into SDP form, whereas the EH constraint is approximated 
by first-order Taylor series expansion. An initialization value 
of the vector u; is given by random generation, and u/ can be 
updated at each iteration, and uj" +1 ' = holds when the 
algorithm converges, in addition, it is guaranteed to converge 
to a locally optimal solution (quite close to the globally 
optimal solution) li25l . Il26l . 

B. Harvested Energy Maximization 

Now, we turn our attention to harvested energy maximiza¬ 
tion problem (Q}, which can be relaxed as 


max t 

W ,t 

s.t. ||w|| 2 < P, 

min|w ff hf| 2 >t, VI, (8a) 

S fc y 0, Vk, (8b) 

where t is a slack variable for the minimum EH constraint, S 
has been defined in 0. Although the secrecy rate constraint 
( l9al ) is LMI, the above problem is still not convex in terms 
of the harvested energy constrai nt d9bl) . The approximation 
is similar to Lemma 1 in Section IIII-A1 in which the energy 

















maximization problem can be reformulated as 

max t 

W ,t 

s.t. ||w|| 2 < P, 

S fc P 0, Vk, (9a) 

||uj”- 1 1| 2 + 2 uj n) (i)[ui(i) — U; n) (*)] > t, VI, (9b) 

i=1 

where u; has been defined in Lemma [T] The above problem 

(n) 

is convex for uj ' at each iteration, and can be solved by 
interior-point method. 


The problem ( ITOt is not convex in terms of channel uncer¬ 
tainties, and cannot be solved directly, thus, we first c onsider 
the reformulation of the secrecy rate constraint (1 1 Oab . which 
can be relaxed as 


i 

<7 S 


w H h s — e s 



> y/ti, 


[ ^(he,fc+e ei fc)' H w ( 2 ^* —1)2] 


^w H (h 

<?e v 

~>R 


kV&e,k) 

(2 k -1)5 


<t 2 , 

(ID 


The first constraint in (ITU is rewritten based on first-order 
Taylor approximation as follows: 


IV. Robust Secrecy Rate Optimization with 
Channel Uncertainty 

In the previous section, we have solve the secrecy rate opti¬ 
mization problem based on the assumption that the transmitter 
has perfect CSI of the eavesdroppers. However, it is not always 
possible to have perfect CSI due to lack of cooperation as 
well as channel estimation and quantization error. Thus, we 
consider robust scheme based on the worst-case secrecy rate 
by incorporating channel uncertainties. 


A. Channel Uncertainty 

In this subsection, it is assumed that the channel state 
information is not available at the legitimate transmitter. The 
channel uncertainties based on worst-case scheme is modelled 
as 


h s — hg-pGs, 

h e,k = he,fc , 

h; = hj+e;, 

where h s , h f /- and h? denote the estimated channels of the 
legitimate user, the fc-th eavesdropper and the Z-th EH receiver, 
and e s , e e j. and c/ represent the corresponding channel errors, 
which are assumed to be bounded as 


—5R{w ff h 3 }- — || w || > /("l(t 2 ), 
(J s (7s 


( 12 ) 


where /( n) (i 2 ) = Jt ( 2 ] H- j===(t 2 - t^). 

2 v4" 5 

Lemma 2: The second constraint in 03 can be reformu¬ 
lated as 


Sk = 


where 


0 

-1 


Sfe.i — \k [ 0 —1 

-E e ,k[dL w H 0 0 ] 


£e,k 


R 

2Z.w H 

0 

0 


Afcl 


yo,Vk. 

(13) 


Sfc.i = 


f (n) (t 2 )I 

R -1 H 

(2 fl — 1)3 


R 

^w H h, 


e,fc 


( 2 r — 1)5 

f (n) (t 2 ) 


(14) 


Proof: Please refer to Appendix [B] ■ 

Thus, the robust power minimization problem can be written 
as 


||e s || 2 = ||h s — h s || 2 <e s ,for e k > 0, 

||Oe,fe||2 — |[he,fe he,fc || 2 ^ , for Ek ^ 0, 

||e(|| 2 = ||h/— h; || 2 < £;, for e; > 0. 

In the following subsections, we will reformulate the secrecy 
rate optimization frameworks based on this channel uncertain¬ 
ties. 

B. Robust Secrecy Rate Optimization with Channel Uncer¬ 
tainty 

Now, we solve the robust power minimization problem by 
incorporating with channel uncertainties, where the CSI of the 
legitimate user, the eavesdroppers and the EH receivers are not 
available at the legitimate transmitter. Thus, the relaxed robust 
problem is formulated as follows: 


min s 2 , 

®2,w ,\ k 

s.t. S2 y K 0, 
w — ’ 

S k (x k ,f^(t 2 ))h0, 

—W ff h s - — ||w|| > f (n \t 2 ), Vk, 

CT s & s 

Klhf 1 w} > E 5 +ei||w|| 2 , w} = 0, VI. (15) 

(n) 

The above problem is convex for a given t 2 at each iteration. 
Thus, an initialization of t 2 is given to solve the problem (1 15 1 ) 
by using interior-point method, which is updated iteratively. It 
is easily observed that t 2 is updated to t < £ l+1 ' > = t 2 when the 
algorithm converges. 


s.t. min log 1 + 


| (h s +e s ) ff w| : 


1 f 1 | (1^6 ,k + ®e,/c) W| \ - j-, w; .... . 

max log 1+ —— 1 --— >R, Vk, (10a) 

e e,k \ at ) 

|(hi + ei) H w\ 2 > E, VI. 


(10b) 


C. Robust Harvested Energy Maximization with Channel Un¬ 
certainty 

Now, we turn our attention to robust harvested energy 
maximization framework, where we consider maximizing the 
harvested energy subject to the achieved secrecy rate and the 
transmit power constraints incorporating channel uncertainties. 






















This optimization framework is formulated as 
max min | (h; + ei) H w| 2 , 

w e; 


s.t. 


min log 

e s 


/ | |(h a + e s ) g w| 2 
V a s 


—max log 1 + 

e e,fc 


|(h ej fc+e ei i ! ) J 


< P, VL 


>R , V&;,(16a) 
(16b) 


The above problem is not convex in terms of channel 
uncertainties and nonconvexity of the secrecy rate constraint, 
and cannot be solved directly. Based on the equivalent refor¬ 
mulations, we propose a two-level SCA based optimization 
algorithm to obtain the robust energy maximization problem 
dl6l) . We denote w* as the optimal solution to the problem 

m. and define — +e ^ —— = r*, then w is also optimal 

solution of the following problem with r = r*. By employing 
a slack variable r, we have 


Fro m Nem irovski lemma shown in Appendix QT] the con¬ 
straint (l2lab can be written as 


At = 


where 


0 

-1 


£e,k 


Ak—fik [ 0 — 1 
-£e,k [ ^-w H 0 o] Mfc 


R 


h 0, (22) 


A*; = 


y/K 

R 1 H 

e , fc 

(2 R — 1)3 


^W"h e ,fc 

e l 

( 2 R -1)2 

y/r 


(23) 


Thus, the robust harvested energy maximization problem is 
equivalently reformulated as 


max min |(h; + ei) H w\ 2 , 

w e; 

s.t. 2 fll(he,fc+e e , fc ) g w| 2 + R 

(T^ 
u e 

|(h s +e s ) ff w| 2 
- 2 -- r ’ 

0-f 

||w|| 2 < P, 


1) < r, Vfe, (17a) 

(17b) 

(17c) 


The proof for statement is easy and the details are omitted 
for brevity. Then, we consider to denote the optimal value of 
the problem (f24ll as /(r), which is a function of r, then the 
following theorem holds, 


Theorem 1: The problem (fT6t is equivalent to the following 
problem: 


max /(t) (18) 

T>0 

Proof: Please refer to Append ix ICl ■ 

Based on Theorem 1, the problem (118b can be solved instead 
of solving ( IT6b . Thus, we propose a two-level optimization 
algorithm to obtain the optimal value of r. In this following, 
we focus on reformulation of the problem (ITTl i based on a 
given r, which can rewritten by employing a slack variable 
for the objective function of (1 1 7b as follows: 


max t 3 , 

W ,t 3 

s.t. |(h; + eij^wl 2 > £ 3 , VI, 

2 n \(h e , k +ee, fc ) w w| 2 + (2 * _ 1} < Tj VJfci 

|(h s + e )"w| 2 ^ T; y 
o’! 

||w|| 2 < P, (19) 

Energy constraint can be approximated based on SCA as 

K{w w h ; } — e;||w ||2 > \jt'£ l> + ■ 


(n) 1 


-.(t 3 -t^),yi. (20) 


2 \ t\ 


The secrecy rate constraint can be reformulated as 

<%/r, Vfc, (21a) 


r r _ 

^w ff (h e ,t + e e ,fc) 


1 


(2 K -1)* 


(K{w H h s } — Es ||w||) > y/r. 


(21b) 


max t 3 

W ,t 3 

s.t. ||w || 2 < P, 

m, I22J, d2lbt . (24) 

The problem (l24l i is convex for a given r and a initialized 
value of the SCA algorithm also can be employed to 
approximate for (12 1 bb . whereas we propose a straightfor¬ 
ward line based search algorithm to obtain optimal r. In order 
to carry out this algorithm, we will determine the iterative 
and feasible region of r (i.e., r m i„ < r < T max ), where 
the optimal value of r* is inside this region. According to 
the problem (ITSl) . we can easily determine the lower bound 
of Tmin = 0- Then, we will show the upper bound r max by 
using inequality properties: 

|(h a +e s ) H w | 2 ||w || 2 112 

-o- < — 5 — h s +e 3 

< —j(||h s ||2 + ||e s H 2) 2 

<?i 

< -^-(||h s || 2 + e s ) 2 . (25) 

The upper bound in ( l25l ) can be derived based on Cauchy- 
schwarz and triangle inequalities, as well as ||w|| 2 < P and 
||e s || < e s . Hence, the iterative region can be obtain as 
0 < t < -^-(||h s ||2 + e s ) 2 - From the reformulation for the 
original problem (fT6b and iterative region, we summarize one¬ 
dimensional search algorithm as Table [TT] 

V. Secure SWIPT of Joint Beamforming and 
Artificial Noise Design 

In the previous section, we solved the secrecy rate optimiza¬ 
tion problems to optimize the transmit beamforming without 
AN, where the proposed algorithms are implemented without 
rank-relaxation at the expense of high computational complex¬ 
ity. In this section, we extend our attention to joint signal 
beamforming and artificial noise (AN) design in secrecy rate 
optimization problem, where the legitimate transmitter sends 
the signal with AN in order to introduce more interferences 
to the eavesdroppers. In addition, AN can be harvested with 
signal beamforming by the EH receivers, which improves 
the performance of SWIPT. The secrecy rate optimization 
problem can be formulated into SDP. Unlike li22l . where it 






























TABLE II: One-Dimensional Search for f24l i 

1) Given Tm ax , At, and an initialization value of 1 = tg°\ 

2 ) Outer Iteration loop begin 

If T < Tmax. then 

a) Inner Iteration loop begin 

If the problem 1241 is feasible, then 

i) Solving G3J based on SCA algorithm for a initialization tg 
based on a given r. 

ii) Set t 3 = tg"). 

else 

Break 

end 

b) Inner Iteration loop end until the SCA algorithm converges. 

c) Updating r = r + At. 

end 

3) Iteration loop end until the required accuracy. 

4) Obtain r* by solving a sequence of problems G! for a given r, and 
the optimal transmit beamforming w* can be obtained. 


1) SDP-Based One-Dimensional Line Search Metho d: We 
first introduce a slack variable t, and the problem (l26l > is 
rewritten as 


max R s +\os(t), 

Qs,w,t 


s.t. 


log 



+°t) 


< log(^), Vfc, 


(27a) 


Tr(Q s + W) < P, Tr[Ai(Q s + W)] < Pi , MR 27b) 
hf(Q s +W)h i>E,,Ml, (27c) 

Q s h o, W >z o,t > 0, 

where A, = is given antenna design parameters to 

adjust each antenna power budget, and a, is unit v'-th vector 
(i.e., [a,;]j = 1 for i = j and [a,]., = 0 for i ^ j ). The specific 
applications of per-antenna power constraint have been already 
described in fl6l . HU. T he p roblem (fTTI i is still not convex 
in terms of the constraint d27ab . and cannot be solved directly, 
thus, it can be formulated as a two-stage optimization problem, 
the outer problem is a function of t, which can be written as 


has shown that the rank of the optimal solution obtained from 
line search algorithm was less than or equal to 2. In this 
paper, we provide a novel SDP relaxation for the secrecy rate 
optimization problem, and exactly show the optimal solution 
always returns rank-one. 


A. Problem Formulation 

We first write the secrecy rate maximization problem with 
the transmit power and the minimum EH constraints, where the 
secured transmit beamforming and AN are jointly designed. 
This optimization problem can be formulated as 


max mini? s — Mk, 

Q.W k 

s.t. Tr(Q s + W) < P, [Q s + W]( i(i ) < Pi, Vi, 

min hf (Q s + W)h, > E u Ml. (26) 

where Q s = E{ww f) } is transmit covariance matrix, W = 
E{vv h } is AN beamforming, [Q s +W]( i i ) ( i = 1 
represents each antenna transmit power constraint, and the mu¬ 
tual information at the legitimate user and fc-th eavesdropper 
can be written as 


R s = log 
Re,k = log 



hfQ s hf 

hFWhF + al 



h^fcQsh 5= \ 
h^Wh l k Pal)- 


Remark 1: The relaxation method to be proposed is also 
suitable for the scenario that the eavesdroppers and the EH 
receivers are equipped with multiple antennas. We will provide 
the analysis in Section [V-DI 


B. Secrecy Rate Optimization 

For the secrecy rate maximization problem d26l >. we will 
provide two different methods to joint optimize transmit 
covariance matrix and AN beamforming: 1) SDP-based one¬ 
dimensional line search, 2) SDP based on SCA. 


max 

t 


log(l + /(f))+log(t), 

S.t. t min ^ t ^ 1, 

The lower bound t m \ u can be derived as follows: 


(28) 


t > 1 + 


h?Q s h J 


> 1 + 


h?Wh« + a\ 

^maxjQsjllhsI 


> 1 + 


hfQ.hf 


> 1 + 


Tr(Q 3 


> 1 + 


P||h 3 


= fn 


(29) 


Then, the inner problem can recast for a given t as follows: 

hfQ s hf 


f(t) = max 
v y q s ,w ,t 

s.t. 


h"Wh» + < 


era - ozD- oo) 

It easily verified that the constraint (1271) can be reformulated 
as 


h H 

n e,fe 


Qs —(- —i)w 


h e , fc <(i-l )*l 


(31) 


Thus the problem can recast for a given t as follows: 

hfQshf 


fit) = max 
Q s ,w,t 

s.t. 


hfWhF + erf ’ 

ob, m, cm . 


(32) 

The problem (f32l> is a quasi-convex problem, thus we consider 
Charnes-Cooper transformation to convert it into a convex 
problem. Let us introduce S so that the following relations 
hold: 


Qs = 


Qs 


W = 


W 


(33) 


Thus, we have 

/(f) = max hfQ s h s , 

Qs,w 

s.t. hfWh s + Sal = 1, 


Qs —( l)W 


h e , k <(±-l)6al 


Tr(Q s + W) < SP, Tr[A,(Q s + W)] < S Pi , Mi, 
h,"(Qs + W)h ; > SE t , Ml, Qs h 0, W h 0. (34) 

The problem ( l34l ) is a convex problem, and can be solved 
efficiently by using interior-point method l27l . Thus, the 





















optimal solution of the problem d32) can be obtained through 
d33l >. once the problem (f34b has been solved. 

2) Tightness Analyses of Rank Relaxation: Now we inves¬ 
tigate the tightness of the rank relaxation to the problem (f32l) . 
we assume that the optimal value f(t ) can be obtained by the 
optimal solution of the problem (I32l> . we have the following 
inequality, 

h«Wh s +af >/(*)^ h f [Qs~f(t)W]h s >f(t)o 2 s , (35) 

Thus, we consider the following power minimization, 
min Tr(Q s ) 

Qs 

s .t. OS), GIB, GIB- 06) 

It is easily verified that the feasible solution of problem (f36l> 
is the optimal of i PH) due to the constraints d35l) . (1 27h t > . (127cl) . 
Thus, the following theorem is provided to show every optimal 
solution of the problem (l36l ) is rank-one: 

Theorem 2: If the problem ( l32l ) is feasible, then there al¬ 
ways exists optimal solution (i.e., Q s ) satisfies rank(Q s ) < 1. 

Proof: Please refer to Appendix [D] ■ 

From Theorem 2, tightness analysis has been provided so that 
the problem (132) has a rank-one solution for all feasible t. 

3) SDP Based Successive Convex Approximation: In this 
section, we propose a SDP based successive linear approx¬ 
imation (SCA) algorithm to joint optimization for transmit 
beamforming and AN, thus, rewrite the problem (f26l) as 

(a e 2 +Tr[h e , fc h5(Q s +W)A L 2 +Tv(h s hfW)\ 

min max -— 

Qs ’ W k (a|+Tr[h s hff(Q s +W)]jUl + Tr(h e , fc h" fc W)J 

s.t. Tr(Q s + W) < P, Tr[A,(Q s + W)] < pi, Vi, (37a) 
hf(Q s +W)h,>£i, VI. (37b) 

The above problem is not convex in terms of the objective 
function. Let us introduce the following exponential variables 
to equivalently modified the objective function. 

e*° = a 2 +Tr[h s hf(Q s +W)], e Xk =a 2 +Tr(h e , fc h" fc W), 
e Vk =<T 2 +Tr[h e ,fch^ fc (Q s +W)], e yo =<7 2 +Tr(h s hf W). 

(38) 

Thus, the problem (l37l > is rewritten by introducing a slack 
variable r as 

min r (39a) 

Qs,XV,x 0 ,y 0 ,x k ,y k 

s.t. e TO ^°" h/fc ^ fc < r, (39b) 

a 2 +Tr[h s hf (Qs+Wjj^e^ 0 , CT 2 +Tr(h eifc h^W) > e Xk , 

(39c) 

a 2 +Tr[h eifc hf fe (Q s + W)] < e Vk , a 2 +Tr(h s hf W) < e yo , 

(39d) 

GIB- GIB- (39e) 

The above problem is not still convex in terms the constraint 
(139db . Thus, Taylor series exp ansio n (i.e., a x +a x lna(j:— x) < 
a x ) is employed to linearise (139db as follows: 

+Tr[h e> fch^(Q s + W)] < e Vk {y k -y k -\-T), (40a) 

a 2 +Tr(h s hfW)<e M (2/o-yo + l), (40b) 


TABLE III: SCA algorithm for the robust secrecy rate 
maximization problem ( 139b . 

1) Initialize (Q s [0].W[0]) so that I39t is feasible, and given k as the 
tolerance factor for stopping criterion. 

2) Iteration loop begin: 

a) Updating ( x 0 [n], x k [n\ , y 0 [n], y k [n]) by (38). 

b) Solving GD with (xo[n],Xk[k\,yo[n],yk[n]) to obtain 
(Q s [n],W[n]). 

3) Iteration loop end until stopping criterion |r(n+ 1) — r(n)| < k. 


Thus, the secrecy rate maximization problem can be relaxed 
as 

min r 

Q s ,W,x 0 ,y 0 ,Xk,yk,T 

s.t. GIB, GIB- GIB. GIB. GB, v/c. (4i) 


From SCA, the approximation with current optima l so¬ 
lutio n can be updated iteratively until the constraints ( 139cb 
and (139db hold with equality, which implies the problem 
(l37l i is optimally solved. This SCA algorithm is outlined 
as Table [III] ■ The optimal solution obtained by the SDP 
based SCA algorithm at n-the iteration is assumed to 
be (Q* (n), W* (n ), (n) , y* 0 (n), x* k (n) , y* k (n), t* (n)), which 
can achieve a stable point when the SDP based SCA algo¬ 
rithm converges, the proof has been shown in |[28l. N ow, we 
consider the tightness analysis for the problem(B9l> due to 
rank relaxation. It is assumed that (Q*, W*) are the optimal 
solution of the problem (l37l > that are obtained by solving the 
problem (ITit with the SDP based SCA algorithm, and the 
correspondin g sl ack v ariabl es (i.e., Xq, j/q, x k , y^, r*) can be 
obtained by (|38I > and ( 139b| i, respectively. Thus, we consider 
the following problem: 

min Tr(Q s ) 

Qs ,w 

s.t. GIB. GIB- 

a 2 +Tr[h s hf ( Q s + W)] > e x °, n 2 +Tr(h e , fc W) > e Xk , 

o 2 e +Tr[h e , k h^ k (Q s + W)]<e Vk , o 2 + Tr(h s hfW)<e y f Vk. 

(42) 

We also assume that the optimal solutions of the above 
problem can be denoted as (Q S ,W), which are the feasible 
solution of the problem (l37l >. Hence, the objective value f 
is obtained by substituting (Q S ,W) into (f3Tb . and we have 
f < t* , which implies (Q. s . W) is at least the same optimal 
solution to (Q*,W*) for ( 137k Thus, provided the problem 
© is feasible for positive secrecy rate, there always exists 
the optimal solution is rank-one, the proof is similar to that 
of Theorem [2] 

C. Robust Secrecy Rate Optimization 

In the previous section, we have solved secrecy rate max¬ 
imization problem based on the assumption that the legit¬ 
imate transmitter has perfect CSI of the legitimate user, 
the eavesdroppers and the EH receivers. However, it is not 
always possible to have perfect CSI at the transmitter due to 
quantization errors and channel estimation. Therefore, robust 
secrecy rate optimization is employed to jointly optimize 
the transmit beamforming and AN by incorporatin g chan nel 
uncertainties, which have been shown in Section TlV-AI In 
addition, robust power constraint per antenna is considered, 
where the Hermitian positive semidefinite (PSD) matrix A, is 













not available at the legitimate transmitter, thus the true PSD 
matrix can be written as 


Ai = A, + A„ || A;||f <£i, Vi, 


(43) 


where A, € H^ T is the estimated Hermitian PSD ma¬ 
trix, A, is estimated error of the matrix A,;, which can be 
modelled as a spherical uncertainty with a norm-bound e*. 
In the following, we consider one dimensional search and 
successive convex optimization methods to solve the robust 
secrecy rate maximization problem by incorporating with the 
channel uncertainties. 

1) Robust SDP Based One Dimension Search: In this sub¬ 
section, robust SDP based one-dimensional search algorithm 
is proposed to jointly optimize transmit beamforming and 
AN by incorporating with channel uncertainties. Thus, the 
two-level optimization framework have been discussed in the 
previous subsection is employed by incorporating with channel 
uncertainties. Since the outer problem does not involv e the 
channel uncertainties that is similar to Section IV-B1I thus 
we focus on the inner problem for a given t, which can be 
reformulated as 


(h s +e s ) ir Q s (h 3 +e s ) 
Q s ,w,t (h s +e s ) ff W(h s +e s )+CT| 


/(f) = max 


s.t. (h e,k+e etk ) Qs-^j-lJW (he^ + e^fc) 

<(i-l)o-e, 

Tr(Q s + W)<P, maxTr[(Ai +Ai)(Q s +W)| <p i; 

Aj 

(in+ei) 11 (Qs+W)(h l +e l )> El, VI, Q 3 >r 0, W>r 0, t > 0. 


(44) 


The above problem is not convex due to channel uncertainties. 
Thus, we consider S-Procedure to solve this robust secrecy rate 
maximization problem, which can be written as 


/(*) = 


(h s + e s ) H Q s (h s + e s ) 


Q s ,w ,t,X etk ,a, (h s + e s ) fl W (h s + e s ) + <r 3 
s.t. Tr(Q s + W) < P, Tr[A i (Qs+W)]-|-£ i ||Qs-|-W|| F <p it 

Vi, (45a) 


Ae.fcl- [Qs - (j - 1 )W] — [Qs —(7 — 1) W]he,fc 
-h" fc [Qs-(i-l)W] c fe 


>r0, Vk, 


a A I+(Q s +W) _ (Qs+W)h( 
hf(Qs+W) hf(Qs+W)h ,-Ei-aie? 


(45b) 

(45c) 


where c k = -h^ k [(Q s ~-l)W]h e}k +{^-l)al-X etk £l. Let us 

introduce a slack variable r to relax the objective function of 
the problem (l45l >. and by exploiting S-Procedure and Charnes- 
Cooper transformation, we have 


/(f) - _ max r, 

Q s ,W,A S ,/j.s ,A e k,oti 


S.t. 


A S I+Q S Qshs 

hfQs h^Qshs — r — A s £s 


' M-w _ _ -whs l. n 

-hfW -hfWhs-(5o-2 + l-/is£s J - ’ 
A e ,*I - [Qs - (i - 1)W] - [Qs - (\ - 1 ) W]h e , fc 
-Kk [Qs— (i—i)w] _ _ c k 

aT+(_Qs+W) _ _ (Q_s+_W)h ; 

hf(Qs+W) hf(Q.+W)hi-^,-aie? _ 
Tr[Ai(Qs + W)[ + £i ||Qs + W|| F < Vi, 

Tr(Qs + W) < 5P, 


>1 0 , 

>r0, 


where c k = -hf t .[Q 8 -( j-l)W]h e , fc +(5(^-l)gg-A e , fc £g. 

By solving the problem d 46 b . we can obtain the optimal value 
/(f)*, which can be written based on channel uncertainties as 


(h 3 + es) w Q s (hs + e 3 ) f(t)* 

(hs + es)«W(hs + e s ) + <71 - 1 ’ ’ 

=► (hs + e s ) W [Qs - /(f)*W](hs + e s ) > /(f)*a 3 2 , (47) 

Thus, the associated power minimization problem can be 
given as follows: 


s.t. 


min 


Qs, W.qj, /3 s ,A eifc 

145ab — l !45cb . 


Tr(Qs), 


/ 3 3 I+[Qs— /(f)*W| [Qs-/(f)*W]hs 
hf[Qs-/(f)*W] d 3 


(48a) 
>: 0, (48b) 


where d s = hf [Q s ~ f(t)*W]h s -~/3 s e 2 s . It easily 
verified that the feasible solution to problem (l48t is optimal 
for (l45l >. which is derived from (I48ab and (I48bb . Thus, the 
following theorem holds to show the optimal solution of fiol i 
is rank-one: 

Theorem 3: Provided that the problem (l45t is feasible, the 
optimal solution of ( l45t always return rank-one, and this 
optimal solution can be obtained by solving (l46t . 

Proof: Please refer to Appendix [E] ■ 

2) Robust SDP Based Successive Convex Optimization: 
Now, we consider the second reformulation for the secrecy 
rate maximization problem to joint optimization for transmit 
beamforming and AN covariance matrix, the optimization 
framework can also be reformulated into robust SDP based 
SCA by incorporating with channel uncertainty. Thus, this 
robust secrecy rate maximization problem can be rewritten as 


min max ^ e ’ fcr ~ s (49a) 

Q s ,w k t s 

s.t. Tr(Q s +W )<P, Tr[(Ai +A;)(Qs+W)] <Pi, Vi, 

(49b) 

hf(Qs + W)h 1 > Ei, VI. (49c) 


where t e ,k = <J^+(h e ^ k +e e ^ k ) H (Qs+W)(h eife +e e ,k}, r s = 
crf+(h s 4e s ) H W(h s -|-e s ), t s = a^+(h s +e s ) H (Q s +W)(h 8 -+e 8 ) 
and r e ,k = cr^+(h. e .k+e e . k ) ^W (h e .k+e e . k ). Let us introduce 
the following relations for ( I49ab 


e x ° < (Js + min(h 3 + e s ) H (Q s + W)(h 3 + e 3 ), (50a) 

e s 

e Xk < er 2 + min(h e! fc + e e ,k) H W(h e , fc + e e ,fc), (50b) 

e e,fc 

e Vk > o-e + max(h e ,fc + e e! i,) H (Qs + W)(h e ,fc + e e ,i,)j(50c) 

e e,fc 

e vo > cr 3 + max(h 3 + e s )' H W(hs + e 3 ), (50d) 

e s 


By emp loyin g the sl ack variables (i.e., r, u s , u e ,fc, v s . and 
v etk ) for (149ab , (I50ab - (l50db . respectively, the problem ( 149 b can 


(46) 
























be equivalently modified as 

min r, 
n 

^yo+yk-^o^k < Tj j49bb . 149c 1 . 


receivers are equipped with multiple antennas. 


s.t. e 
x 0 J 2 


e°°° <a s +u a , min(h s +e s ) [Q s +W](h s +e s ) >m s , (51a) 

e s 

e Xk <Ve+u e ,k, min(h ej fe+e ej fc) ir W(h Ej fc +e e ^)> u et k, Vk, 

e e,fc 

(51b) 

e 2 "" >cr 2 +v e ,k, max(h e ,fe4-e e ^)' H (Q s + W)(h e ,fe +e e ^)<D e .fc, 

e e,fc 

(51c) 

e M >oy+^s, max(h s + e s ) ff W(h s +e s ) < v a , (5Id) 

e s 

{Qsi W, e s , e e ,fc ,Xo,yo,Xk,yk,U s ,Ue,k,Vs,Ve : k} £ 9. (51e) 


By exploiting S-Procedure and first-order Taylor series 
approximation, we have 


mm t, 
n 


bO, 


s.t. e voiyk ~^ 0 ^ k <r, e x ° < a 2 +u s , e Xk <al+u e} k, 
e Vk (yk-yk + l)>o-e+v e ,k, e yo (yo-yo + l)>Vs+v s , 

~ A £ I+(Q S +W) _ (Q s +W)h s 

h?(Q s +W) h?(Q.+W)h,-u.-A,e2 
Ae,fcl + W _ _ Wh e ,fc 

h^W h^ k Wh etk -u e , k -Xe,ksl,k _ 
3 e ,fcI-(Q s +W) _ -(Q s +W)h e , fc 
-h" fc (Q s +W) -h^ fc (Q s +W)h e,k+V ei k-Pe,k£ 2 e,k 


t o, 


S I w 


-Wh s 


-hfw -h»Wh s +v s -l3s£ 2 s 


(52a) 

(52b) 

(52c) 
(5 2d) 

to, 

(52e) 

(52f) 


log 

Tr 


1+ a'I+H" t WH e , t H" fe Q s H e , 


Hf(Q s +W)H ; 


> Ei. 


<log (t 1 ), (53a) 


(53b) 


where H € C NTXNE ’ k and H; G qN t xNi ^g channel 
matrices between the transmitter and the fc- th ea vesdropper as 
well as the Z-th EH receiver. The constraint (I53ab can be easily 
relaxed as 


{\ ~ l)a e 2 I + H lk{{\ ~ 1)W - Q s ]H e , fc X 0. 


(54) 


QiI+(Q s +W) 


(Q s +W)h( 


Tr(Q s +W )<P, Tr[A4Q s +W)] + ||Q s +W || F < Pi , 

[Qs X 0, W X 0, X 0 , yo, Xk , yk , Us , Ue,k ,V a ,V B ,k, 

A s > 0, A e ,k > 0, /3s > 0, Pe,k > 0, at > 0] G Q, Vi, l, k. 


(52h) 


The above problem is convex for a given y k and y () at 
each iteration, can be solved by using interior-point method to 
update the solution for the next iteration until the algorithm 
converges. Thus the robust SDP based SCA algorithm is 
similar to Table [HI] On the other hand, the rank -1 solution 
can be obtained that is similar to that of one-dimensional line 
search algorithm shown in Theorem [3] 


This proof for the above LMI is omitted due to space limit, 
readers are referred to m for more details. Since (l54l i is 
obtained based on the as sump tion that rank(Q s ) < 1, we 
substitute t he co nstra ints (I53bb and (l54l > into the associated 
constraints (127al l and ( 127cb in the Lagrange dual function ( l68l >. 
where it easily verified that the optimal solution of Q s always 
returns rank-one. Hence, the problem ( l32l > still has rank-one 
solution when the eavesdroppers and the EH receivers consist 
of multiple antennas. Secondly, we provide a proof for the 
optimal solution of il45l i still returns rank-one for multiantenna 
eavesdroppers and EH receivers. The channel uncertainties can 
be modelled as 

H ei fc = H ei fc-t-E ei k, Vk, H( = Hi+E ! , VI, 

where H e ^ and H/ are estimated channel, and E e> k and E; 
are channel error matrices, which follow ||E e fe||i? < e k and 
||E;||^ < £;. Thus, the constraints (l53l > can be relaxed by 
incorporating with channel uncertainties as 


(55) 


where Ai = 

A 2 = H^ fc ((f- 1 -l)W-Q s ), A 3 =\t~ 1 - 1)W-Q S +^I, 

B! = a z I+[I®(Q s +W)], B 2 = [I® (Q a +W)]hj, and B 3 = 
hp[l 8 )(Q s -(-W)]h;— Ei—aief. The first LMI of ( 133 b is obtained 
based on M Lemma 5], whereas the second LMI can be 
relaxed by exploiting S-Procedure. Thus, these two constraints 
are replaced the corresponding constraints in the dual function 
(ED so that rank-one solution still can be obtained for the 
robust secrecy rate maximization problem (l45l >. ■ 

VI. Numerical Results 





Ai 

A 2 

^ 0 , 

Bi 
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aiSi 

to, 

(52g) 

. A ? 

A 3 

. B 2 

B 3 


D. Extension to The Case of Multiantenna Eavesdroppers and 
EH Receivers 

Now, we show that the Theorem [2] and Theorem [3] can 
also be applied in a more challenging scenario that the 
eavesdroppers and EH receiver are equipped with multiple 
antennas by the following corollary. 

Corollary 1: If the eavesdroppers and EH receivers are 
equipped with multiple antennas, the optimal solution of the 
problems d32l ) and < 143b still always return rank-one based on 
perfect and imperfect CSI. 

Proof: We first reconsider the constraints (127ab and < 1 27c b 
based on the assumption that the eavesdroppers and EH 


In this section, we provide the simulation results to validate 
our proposed algorithm. We consider the a MISO secrecy 
system in the presence of three eavesdroppers and two EH 
receivers. The legitimate transmitter is equipped with four 
transmit antennas (i.e.. Nr = 4), whereas the others are 
equipped with single antenna. All of channel coefficients (i.e., 
h s , h e y, and h/) have been generated as zero-mean circularly 
symmetric independent and identically distributed Gaussian 
random variables. The noise power are assumed to be 1 (i.e., 
of = <Tg = 1). The transmit power (i.e., P) is assumed to 
be 20dB and all of the error bounds (i.e., e s , e e and £/) are 
replaced by e for convenience, which is set to be 0.05 unless 
specified. 



















A. Secrecy Rate Optimization without AN 

We first provide the simulation results for the secrecy rate 
maximization problem for perfect and imperfect CSI. Fig. Q] 
shows that the achieved secrecy rate with different transmit 
powers for robust schemes, non-robust schemes and perfect 
CSI case, where it is easily observed that the achieved secrecy 
rate increases with transmit power, and our proposed scheme 
achieves the same performance with the scheme with SDP via 
rank-relaxation. In addition, the robust scheme outperforms 
the non-robust scheme in terms of achieved secrecy rate. 



Fig. 1: Secrecy rate with different transmit powers. 


Then, the harvested power with different transmit powers and 
different error bound e are shown in Fig. [2] and Fig. 0 From 
these two results, the harvested power increases with the 
transmit power and decrease with the error bound (i.e., e), 
respectively. 



Fig. 2: Harvested energy with different powers 



Fig. 3: Harvested energy with different error bound e 


B. Secrecy Rate Optimization with AN 

In this subsection, we provide the simulation results to 
validate the proposed secrecy rate maximization problem for 
joint optimization of transmit beamforming with AN. Fig. 
Q] shows that the achieved secrecy rate with different the 
ratio of the harvested power with the transmit power (i.e., 
77 = based on perfect and imperfect CSI, where one¬ 
dimensional line search based scheme and SCA algorithm with 
or without AN are plotted, respectively. From this figure, the 
achieved secrecy rate increase with transmit power, and both 
schemes in perfect CSI has the same secrecy performance 
with or without AN, whereas the SCA based scheme has 
a little better performance than one-dimensional line search 
based scheme at 0-25 dB, the SCA based scheme outperforms 
slightly that with one-dimensional line search from 25 to 
30 dB. Moreover, the transmit beamforming with AN has 
a better performance than that without AN in terms of the 
achieved secrecy rate. Next, the convergence performance of 
the SCA based scheme is discussed in Fig. [5] which shows 
the convergence of the SCA scheme, from this figure, we 
can observe that the SCA algorithm converges slower as the 
transmit power increases. Then, the achieved secrecy rate with 
the ratio of the harvested power E to the transmit power 
P is presented in Fig. [ 6 ] where we compare our proposed 
schemes (i.e., one-dimensional line search and SCA) with 
two-dimensional search and the scheme shown in l22l . It is 
observed that one-dimensional line search scheme is closing 
to the scheme in ll22ll and two-dimensional search scheme. 
In addition, the SCA based scheme outperforms these three 
schemes. At last, we provide the simulation result to show 
the proportion of the AN power Tr(W) in the total transmit 
power ( P ) with different signal-to-noise ratio (SNR), which 
can be replaced by P due to the noise power equals to 1. One 
can observe that this proportion increase first and then has 
a decline after approximately 15 dB, and the robust scheme 
without the EH receivers has a lower proportion of the AN 
power than the robust scheme with the EH receivers at the 













































lower SNR regime. In addition, the scheme in lf22ll has a larger 
proportion than our both proposed schemes. 



Fig. 4: The achieved secrecy rate with different transmit power 



Fig. 5: The convergence performance of SCA algorithm in different 
transmit powers 



Fig. 6: Robust secrecy rate with different the ratio of harvested 
power with transmit power |j. 



Fig. 7: The proportion of AN consumption with different SNRs. 

problem exactly returns rank-one. Simulation results have been 
provided to validate the performance of our proposed schemes. 


VII. Conclusions 

In this paper, we investigated the secrecy rate optimization 
problems for a MISO secrecy channel for SWIPT in presence 
of multiple eavesdroppers and EH receivers. We first proposed 
a novel reformulation by using Nemiroski lemma to optimize 
the transmit beamforming to circumvent rank-relaxation based 
on perfect CSI and imprefect CSI. Then, jointly optimization 
for transmit beamforming with AN are considered in the 
secrecy rate maximization problem based on both perfect 
CSI and imperfect CSI. We reformulated this optimization 
problem by using one-dimensional line search method and 
SCA algorithm, respectively. Furthermore, we provide the 
tightness analysis for this optimization problem, which can 
show the optimal solution of this secrecy rate maximization 


Appendix 


A. Proof of Lemma [7] 

In order to prove Lemma 1, we first rewrite the secrecy rate 
constraint in <[6| as 
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Then, the following lemma is required to convert the above 
constraint into a SDP: 


Lemma 3: (Schur complement) EH: Let X is complex 
hermitian matrix, 
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Thus, S = C — B ff A 1 B is the Schur complement of A in 
X, and the following statements hold: 

• X >~ 0, if and only if A >- 0 and S >*- 0. 

• if A 8 - 0 then X >- 0 if and only if S >- 0. 

By exploiting the Schur complement, (l56b can be reformulated 
as 
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In addition, we consider the reformulation of the EH con¬ 
straint in ©. In order to express this constraint clearly, we 
introduce two variables (i.e., xi £ R. and yi e M) such that 
this constraint can be equivalently modified as 


xf + yf > E, (59a) 

xi = K{w ff hi}, yi =Sr{w H h;}, VI. (59b) 

The constraint d59bl> is convex (linear), whereas ( 159ab is 
not still convex, thus, first-order Taylor approximation is 
considered to obtain the desired upper bound. 

Setting u/ = [a;* yi\ T , thus, xf + if f = u^u;. uf' 1 is n-th 
iteration of the vector u/. Thus, (159ab can be approximated as 
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i= 1 

where i denotes the /-th element of the vector u/. 

This completes Lemma [T] ■ 
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Thus, the robust secrecy rate constraint can be reformulated 
as 
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This completes Lemma [2] 


C. Proof of Theorem [7] 

In order to prove Theorem 1, it is assumed that ipi and ip 2 
are optimal values of ( IT 6 b and ( 1 1 8b , res pec tively. We first show 
(ITSb can obtain the optimal value of d 16b (i.e., tp 2 > T\)■ It 
easily verified that the following equality holds: 

ipi = f(r*). (67) 

On the other hand, ip 2 = max T >o f(r) > f(r*), thus, ip 2 > 
<p\. Secondly, in order to show that ( IT 6 b achieve the optimal 
value of (IT 8 l) (i.e., ipi > ip 2 )- Here, we assume that is the 
optimal value of (fl 8 l) . and w' is the optimal solution of (ITTb 
with t = rf It can be easily observed from d 1 6ab . d 1 7ab . and 
d 1 7bb that wt is a feasible solution of (IT 6 l) . thus, ipi > ip 2 . 
We combine these two parts, it includes tp\ = p 2 . 

This completes Theorem Q] ■ 


B. Proof of Lemma \7} 

The second constraint in (ITTb can be rewritten as 


D. Proof of Theorem \2\ 

The Lagrange function of ( l36b can be written as 
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Lemma 4: (Nemirovski lemma) (29|: For a given set of 
matrices A = A H , B and C, the following linear matrix 
inequality is satisfied: 

A V BXC + C H X H B, ||X|| < t, (62) 

if and only if there exist non-negative real numbers a such 
that 
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By exploiting Nemirovski lemma, the constraint in (16Tb is 
written as 
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where Y, Z, A, p, , r/,, vj , r k denote the dual variables of 
q s , w, m, (127bb . (127cb . and (131b , respectively. Then, we 
consider the following related KKT conditions: 


dC 

dq t 


dC 

dW 


Q S Y 


n t l 

0, =► Y = I-Ah s hf +pl+Y^ ViAi-^2 uih t hf 

i=l 1=1 

K 

+E Tfehe ’ fch ? fe5 (69a) 

k =1 

N t L 

0 , => Z = A/(t)h s hf+/iI+^ ViAi-Y'Vthih? 

i=1 1=1 

K 

-YTk(t-l)K,k<k, (69b) 

k =1 

0,z y 0, A > 0, (69c) 











































By subtracting (I69bb from ( I69a| i. we have 


Y - Z - I — A(1 + /(t))h s hf + r k th e , k h* k , 

k= 1 

=> Y = A —A(l + /(f))h s hf, 


(70) 


By subtracting d72bb from (|72a[ >, we have 
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where A = I + Z + ffi'f 1 'T; 1 ./h,, /.h ( /, ;; .. From (ITOt . one can 

easily observe that A is positive definite, and rank(A) = Nt, u - , ra , ,, 

n wvt at at i tt n n We premultiply d73j by Q s , 

whereas rank(Y) = Nt or Nt — 1. However, Q s = 0 when 

rank(Y) = Nt due to the first condition in (17 2c b . which 

violates R s — R e ^ > 0. In addition, it is easily verified that 

A > 0 and f(t) > 0. Thus, rank(Y) = Nt — 1 always holds, 

which implies Q s lies in the null space of Y from (l72cl i. thus 

rank(Q s ) = 1. 

This completes Theorem 2. ■ 
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The following rank relation holds: 


rank(Qs) = rank 


E. Proof of Theorem [T] 

We first consider the dual function of the problem 
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£(Q S ,W, Y, Z, A,7i,Ts,T ejfc ,T ; )=Tr(Qs)—Tr(YQs) Based on the above rank relation, we need to show 

n t . rank(H s T s H^) < 1 if we claim ran^Q^) < 1, thus, we 
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where Y <E Z G H^ T , A G R+, 7i G R+, T s G H^ T , 

Te.fc £ H? T and T; G are dual variables of Q s , W, 

iaji, (14. Sel l and < 145hb . respectively, in addition. 
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Lemma 5: If a block hermitian matrix P = 


The related KKT conditions are considered as follows: 
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then the main diagonal matrices Pi and P 4 are always 
matrices 11301 . 

Now, we can claim /3 S I+[Q S — /(f) W] ^0 and is nonsingular, 
thus pre(post)multiplying by a nonsingular matrix will not 
change the matrix rank. Thus, the following rank relation 
holds: 

rank(HsTsHf) = rank^ [ Ojv t h s ] T S H 

< rank^ [ 0jv t h s ] ^ < 1. (77) 

This completes Theorem [3] ■ 
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